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. 1 Introduction 

o: 

\ Under a curvature condition, Wang |21j introduced dimensional-free Harnack inequality for 

diffusions on Riemannian manifold. This type of inequality has been studied extensively, see, 
for example, Aida and Kawabi [21 [3] for infinite dimensional diffusion processes; Wang [23] for 
^ \ stochastic generalized porous media equations; Rockner and Wang [17] for generalizes Mehler 

■ semigroup; [1] for stochastic functional differential equation; Ouyang [20] for Ornstein-Uhnelbeck 
processes and multivalued stochastic evolution equations etc. 

Harnack inequality has various applications, see, for instance, [H [T71 [181 1221 (23] for strong 
Feller property and contractivity properties; [21 [3] for short times behaviors of infinite dimen- 
sional diffusions; [HI [11] for heat kernel estimates and entropy-cost inequalities. [21 [121 fl?! [2T] 
established Harnack inequalities using the method of derivative formula. In order to obtain 
Harnack inequality on manifolds with unbounded below curvatures, [5] introduced the approach 
of coupling and Girsanov transformations. In the paper, we will use the above two methods to 
establish Harnack inequalities for stochastic differential equations driven by fractional Brownian 
motion. 

Solutions of the stochastic differential equations driven by fractional Brownian motion have 
been studied intensively in recent years, for example see [131 US] using the pathwise approach; 
see [9] using the tools of rough path analysis introduced in [13]. We prove Harnack inequality 
for stochastic differential equations driven by fractional Brownian motion with Hurst parameter 
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H < ^. As applications of the Harnack inequality, the strong Feller property and the log- 
Harnack inequality are derived. We also get the entropy-cost inequality with respect to the 
Euclidian distance. 

The paper is organized as follows. In section 2, we give some preliminaries on fractional 
Brownian motion. Section 3 prove the Harnack inequality by using the approach of coupling 
and Girsanov transformations, and present their applications. In section 4, we are devoted to 
establish derivative formula and give the corresponding Harnack inequality. 

Harnack inequality 

2 Preliminaries 

Let = {B^^ ,t G [0, T]} be a fractional Brownian motion with Hurst parameter H € (0, 1) 
defined on the probability space Q,, P, i.e., B^ is a centered Gauss process with the covariance 
function 

Rnit, s) = E(Sf ) = i {t'^ + s'^ -\t- sD . 

In particular, if H = ^, B is a Brownian motion. It is well known that if H ^ ^, B^ does not 
have independent increments and has a-order Holder continuous path for all a £ {0,H). 

For each t € [0, T], we denote by Tt the cj-algebra generated by the random variables {B^ : 
s G [0,T]} and the P-null sets. 

We denote by the set of step functions on [0, T]. Let V. be the Hilbert space defined as the 
closure of S' with respect to the scalar product 

{Ilo,t]J[o,s])RH{t,s) = E(i?f ) = 1 {t'^ + s'^ - \t - s\^^) . 

The mapping IjO, T] i-> B^^ can be extended to an isometry between V. and the Gauss space Tii 
associated with B^ . Denote this isometry by i— > B^ {(/)). For more details, one can see |16j . 
On the other hand, from [1^, we know the covariance kernel Rnit, s) can be written as 

rtAs 

RHit,s)= KHit,r)KHis,r)dr, 
Jo 

where Kh is a square integrable kernel given by 

Knit, s) = r{H + ]-r\t - ef~'^FiH H,H + - i), 
2 2 2 2 s 

in which F {■,-,■,■) is the Gauss hypergeometric function. 
Define the linear operator K'^S' — > L^[0, T] as follows 

iK*Hcl>){s) = KHit,s)Hs) + j\ct>{r) - ct>{s))^{r,s)dr. 

By [U, we know that, for all </>, V' G (i^|/-<A, K*jjip)\[{), T] = (0, V') holds. From B.L.T. theorem, 
K"^ can be extended to an isometry between H and L^[0, T]. Therefore, according to [5, the 
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process {Wt = B{{K'^) ^(/[oi])),t G [0, T]} is a Wiener process, and has the following 
integral representation 

Sf = f KH{t,s)dW,. 
Jo 

By [10], the operator Kh : L^[0, T] — > Iq^ ^(L^[0, T]) associated with the square integrable 
kernel Kh{-, •) is defined as follows 

(KHfm := [' KH{t,s)f{s)ds, /GL2[0,r]. 
Jo 

H+- 

where /q^ ^ is the a-order left fractional Riemannian-Liouville integral operator on 0, T, one 
can see [19]. It is an isomorphism and for each / G L^[0,T], 

{KHf){s) = Il_,s"-h^;'^s^~^f, H>\. 

H+- 

As a consequence, for every h G /o_|_ ^(L^[0, T]), the inverse operator Kjj is of the following 
form 

H—- --H 

where I?q_,_ ^ {Dq_^ ) \s H — ^(J^ — f/')-order left-sided Riemannian-Liouville derivative, one also 
can see [T9] . 

In particular, if h is absolutely continuous, we have 

{Kj,'h){s) = s^-hl;''s'^~^h\ H<\. 

In [H], D.Nualart and Y.Ouknine discussed the following stochastic differential equations driven 
by fractional Brownian motion on M, 

dXt = b{t,Xt)dt + dBf^, Xq = X. (2.1) 

They proved the existence and uniqueness of a strong solution for ()2.ip when x) is a Borel 
function with linear growth in x in case H < ^. 

The aim of the paper is to consider the Harnack inequality for the equation (j2.ip in case H < ^. 
We define Ptf{x) := E/(Xf), t G [0,T], / G =^;,(M), where Xf is the solution to the equation 
()2.ip and ^^(M) denotes the set of all bounded measurable functions on M. 

3 Main results and proofs 

Let us start with the following hypothesis (HI): 

(i) \b{t, x) — b{t, y)\ < K\x — y\, Vx, G M, t G [0, T], where -ftT > is a constant; 
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(ii) The mapping t i-^ b{t,0) is bounded on [0,T]. 

It is clear that under (HI), the equation (]2.ip has a unique solution. Furthermore, we can give 
the Harnack inequality for the equation (|2.ip as follows. 



Theorem 3.1 If (HI) holds, then for any nonnegative f € ^{,(M) and t > 0,x,y G 



where C{T, K, H) 



iPrfiyW < PTnx)exp[^CiT,K,H)\x - y\% 

p — 1 

_B(|-£^|-H)\^ r2-2g 



Proof. The proof will be divided into three steps. 
Step 1. Consider the following coupled stochastic differential equation 



dYt =b{t,Yt)dt + dBf + utdt, Yo = y, (3.1) 
dXt = b{t, Xt)dt + dSf , Xo = x, (3.2) 



where the drift term ut of the equation (13. ip is of the following form 

Xt-Yt , 
^*-^^^{*<-}' 

r is the coupling time of Xt and Yt defined by 

r = inf{t > : = Yt}, 

and ?7j is a deterministic function on [0, oo) specified later such that the force ut can make the 
two processes X and Y move together before time T. 

It is obvious that the assumption (HI) implies \h{t,x)\ < C{1 + |x|), then, according to [HI 
Theorem 8], the equation (13. ip has a unique solution. 

Note that d{Xt—Yt) = {b{t, Xt) — b{t, Yt)) dt—utdt, thus applying the Tanaka formula to 1X^—1^1, 
we have for t < t 

d\Xt-Yt\ = sgn{Xt-Yt)d{Xt-Yt) 

= sgn{Xt - Yt) {b{t, Xt) - b{t, Yt)) dt - rjtdt. 

By (HI), for ah t < r we get 

d\Xt-Yt\ < {K\Xt-Yt\-i]t)dt. 

This implies that 

g-K(TAr)|^^ _Y^\<\x-y\+ e-^'r]tdt. (3.3) 

Jo 

Choosing 

Vt = ^ ■ \x - y\, t>0. 

We conclude that t <T and Xt = Yt, a.s. Otherwise, if r > T, by (|3.3p we get Xt = Yt- But 
this contradicts with the assumption that t > T. 
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Step 2. Let B^^ = J^Ugds + B^^ , Vt G [0, T]. By simple calculus, we know that ufdt < 



oo. 



Hence, J^Urdr E /^^^ (L^([0, T])). According to integral representation of fractional Brownian 
motion and the definition of the operator Kh, we deduce 



B 



H 



Usds+ / KH{t,s)dWs 



KH{t,s) 



{Kjj^ / Urdr){s)ds + dWs 



KH{t,s)dWs. 



Now, let 



Rt = exp 



K 



H 



Urdr ] {s)dWs - - 



Urdr 1 {s)ds 



Next we want to show (^/^)o<t<T is an ^j^^ -fractional Brownian motion with Hurst parameter 
H under the new probability RtP- Due to [14^ Theorem 2], it suffices to show that "KRt = 1. 
Since /q Ur-dr is absolutely continuous, then 



S2 U.. 



Hence, we have 



{Kjj^ [ urdr){t 
Jo 



ra-H)' 



2 / r2 ^Ur{s — r) 



2 dr 



< 



< 



ra-H) 



|?7r.|r2 ^(^s — r) ^ 2dr 





\x - y\ 



ra-H) (2K)-i(l-e-2^^ 



■2 [s — r) 2dr 



Ba-H,\-H) 



\x - y\ 



ra-H) {2K)-^{l-e-^^^ 



As a consequence, we get 



Eexp 



K 



— 1 

H 



Urdr I {s)ds 



<exip[C{T,K,H)\x-y\% 



(3.4) 



where C{T, K, H) = (^—^Y{ri^{y~^^ j^,2(i„J-2Kr)2(i_jj--) • Using the Novikov criterion, we have 
ERt = 1. 

Step 3. From step 2, we can rewrite (j3.ip in the following form 

dYt = b{t,Yt)dt + dBl^, Yo = y, 

where {B^)Q<t<T is an ^^^^ -fractional Brownian motion with Hurst parameter H under the 
new probability RtP- By the uniqueness of the solution and Xt = YT,a-s., we have 



Prfiy) = Efixy,) = ERTf{Y^) = ERrfiX^). 



(3.5) 
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Applying the Holder inequality to (I3.5p . we obtain 

iPTf{y)r < PrFix) ■ {ERf~'r~\ (3.6) 
Now we will estimate moments of Rt- 

Denote a = ^ and My = - Jq {KJj^ /g Urdr) {s)dWs. Since {Rt)o<t<T is a P martingale, by 
(|3.4p we have 

ER^ = Eexp[aMT - ^a{M)T] 

= Eexp[aMT - ^a^{M)T + ^a{a - 1){M)t] 

< e^p[a{a-l)C{T,K,H)\x -yl"^]. (3.7) 
Substituting ()3.7p into (|3.6p . we get the desired result. 

Remark 3.2 In the proof of Theorem \3.1\ the choice of ut is not unique. For instance, we can 
take another as follows 

Correspondingly, the result of Theorem I3.il is of the following form 

{PTf{y)Y < PTnx)exp[^C{T,K,H)\x - y\% 

p — 1 

where C{T, K, H) = (^^^^^^^^^ 4K-Hi-e-'^"ni-H) - 

As applications of Theorem 13. H we prove the following results on strong Feller property for 
Pt and log-Harnack inequality. 

Proposition 3.3 Assume (HI). Then Pt is strong feller and the following estimate holds 

\PTf{x)-PTf{y)\ < \\f\\oo[2C{T,K,H)]'2\x - y\exp[C{T,K,H)\x - yW 

for every T > 0,x,y £R and f e ^b(M). 

Proof. It follows from the proof of Theorem 13.11 that . for each / G 

iPrfix) - PTf{y)\ = |E/(X|) - ERTf{X^)\ < ||/||ooIE|l - Rt\. (3.8) 

Next we will estimate the term E|l — Rt\. 
Firstly, we have 

{E\l- Rrlf <E\1- Rt\^ = ER'^-1. (3.9) 

Taking a = 2 in (j3.7p . we have 

ERI <exp[2C{T,K,H)\x -yl"^]. (3.10) 

Combining (j3.9p with (j3.10p . we get 

(E|l - Rrlf < 2C{T,K,H)\x-y\^(i-x.p[2C{T,K,H)\x-y\\ (3.11) 

where we use the elementary inequality — 1 < xe^ , Vx > 0. 
Substituting (13. lip into (j3.8p . we can deduce the desired result. 



Corollary 3.4 Let (HI) hold, then 

PT(log/)(x) < logPrfiy) + C{T,K,H)\x - y\^, 
Vx,y G M,t > 0,/ > 1,/ € =:^b(M). 

That is, log-Harnack inequality holds. 

In fact, since R is a length space, then, by [21^ Proposition 2.2], we know the result holds. 

To state further application of Theorem 13.11 let us introduce another assumption and some 
notations. 

(H2): let /i be a probability measure on R such that for some K > 0, 

MPt/) </W),V/G^+(M). 
Note that if fi is Pr-invariant, then (H2) holds. 

Remark 3.5 The measures jj, satisfying (H2) always exist. For instance, 

oo ^ 
n=l 

where {P^{x, •))n>i is defined recursively as follows 

Pt{x,A):=PtIa{x), P^{x,A):= [ P^-\x,dy)PT{y, A), n>2. 

JR 

Let ^(//, 1^) denote the set of all couplings of /i and u, where /.i and v are two given probability 
on R, and W2{lJ,, y) be the L^-Wasserstein distance between them with respect to the Euclidian 
distance, i.e. 

W2{l^,v)= inf / / \x -y\'^'K{dx,dy). 

Corollary 3.6 Assume that (HI) holds and fj, satisfies (H2) {K = 1). Then the following 
entropy-cost inequality holds for each T > and f € =^^(M) with fJ,{f) = 1, 

^,{P^f\ogP^f) < C{T,K,H)Wiifi,ffi), 

where P^ is the adjoint operator of Pt in L'^{fi). 

Proof. By Corollary 13.41 for P-^f, we have 

Prilog P^f){x)< log PT{P^f){y) + C{T,K,H)\x-y\\ Vx,2/GM. (3.12) 
Integrating both sides of (j3.12p with respect to vr € '^(/i, //x), we get 

li{P^ f log P} f) < lx{\ogPT{P^f)) + C{T,K,H) f [ \x-y\\{dx,dy). 

jR JR 

Note that, the Jensen inequality and the hypotheses imply 

fi{logPT{P^f)) < log fi{PT{P}f)) < log fi{P^f) = logfiifPrl) = log/i(/) = 0. 
So, we get 

fi{P^flogP^f)<C{T,K,H) inf / / \x-y\\{dx,dy). 

7re^(M./W JR JR 

The proof is complete. 
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4 Derivative formula 



In this part, we begin with the following hypothesis (H3): 

(i) d2b{t,x) < if, Vx € M, t G [0, T], where K > is a constant, where d2b{t,x) denotes the 
derivative for the second variable; 

(ii) The mapping 1 1— )• b{t, 0) is bounded on [0, T]. 



The aim is to establish a Bismut type derivative formula for Pt which will imply the Harnack 
inequality. For / G i^;,(M), x, y G M, T > 0, we will consider 



DyPrfix) := lim 



Prfix + ey) - Prfix) 



Theorem 4.1 (Derivative formula) Assume (H3). Then, for each T > 0, / G ■^^b(]R),x,y G 
R, DyPTf{x) exists and satisfies 

DyPrfix) = Ef{X^)NT, 



where Nj- 



r{i-H)T io 



lo 7^^(1 + d2b{r,Xr){T - T))dr 



(s-r-)3 



ydWs 



Proof. As above, solves the equation (j2.ip . For any e > and y G M, we introduce the 
following equation 



dXl = b{t, Xt)dt + dB^ - ly, = X + 



(4.1) 



By (H3), we easily know that the above equation has a unique solution. Combining (12. Ih 
with (j4.ip . we deduce that Xf — Xt = ^^ey, Vt G [0, T], in particular, X^ = Xt- Let 
r]t = b{t,Xt) - b{t,Xl) - j^y, G [0,r], then we can rewrite (f4T]l in the form: 

dXl = b{t,Xt)dt + dBf, 



where = Bp + r/^ds. Note that 

— , ^, ey — T — t ev 

\Vt\ < K\Xt - Xl\ + ^ = K^ey + ^, 

SO, we have /g rjjdt < oo, and moreover, j^rfrdr G /q_,_ ^ (L^([0, T])). Due to the integral 
representation of fractional Brownian motion and the definition of the operator Kh, we get 



KH{t,s)dWs 



where Wt = Wt + jl{KJj^ J^r]rdr){s)ds. Now, let 



Re = exp 



K 



-1 

H 



i]rdr I {s)dWs 



T 



K^j/ / Vrdr j {s)ds 
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Now we will prove that (i?^ )o<t<T is an -fractional Brownian motion with Hurst parameter 
H under the new probability ReP, according to \\.4\ Theorem 2], it only needs to show Ei?^ = 1. 
Similar to step 2 of theorem 13.11 we get 



< 



3 _ tr 1 
2 -"'2 



ra-H) t' 



Hence, it follows that 



Eexp 



T 



K^^ J rjrdr j {s)ds 



< oo. 



By the Novikov criterion, MRt = 1 holds. 

Hence, in view of the uniqueness of the solution and = Xt, we have 

PTf{x + ey)=ERJ{Xf). 
With the help of the dominated convergence theorem due to ()4.2p . we deduce that 

Prfix + ey) - Prfix) 



DyPrfix) 



lim 

e-i-0 

lim 

e->0 



E 



e 

Re - r 
I 

e 

R,-l 



/(X|)lim 



(4.2) 



Let Mt =: - /o^ [KJj^ /q r/^dr) {s)dWs. Thanks to we get 

{M)t= r \{Kh^ [ iirdr){s)\^ds<Ce\ 
Jo Jo 

where C is a positive constant. Therefore, we deduce 



Re-1 

hm 

e-^0 € 



exp[MT - ^{M)t] - 1 
Mt-\{M)t 



lim 

lim 

Mt 
iim . 

e-fO e 



Note that 







1 



Mt = - I {K^' 1^ r]rdrj {s)dWs 

,-, . s 2 T-rirdrdWs 

r(i - H) Jo Jo {s - 

s^"^ / ^[h(r,Xl) -h{r,Xr]drdWs 



1 



va-H) 







(s-r)^+2 



+ e 



n\-H) Jo 



T 



T 2 



-H 



y 



(s — r) "^2 



T^drdWs, 
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therefore by (H3), we conclude that 



hm 



R,-l 



1 



T 



e Ti^-H)Tjo 
The proof is complete. 



1_ rr 
f2 



(s — r)2 



{l + d2b{r,Xr){T -r))dr 



ydWs. 



Remark 4.2 IfH=^, i.e. is a Brownian motion, then the corresponding derivative formula 
is of the following type 



DyPrfix) = i /^[l + (T - s)d2b{s,Xs)]ydWs 
Jo 



Remark 4.3 Since we deal with one dimensional case, the derivative formula of theorem ()4.ip 
is equivalent to Prfi') is derivative. The method we adopt is also valid for n-dimensional case. 

As an application of the derivative formula derived above, we have the following result. 



Corollary 4.4 If (H3) holds, then for any nonnegative f € 

P 



and t > 0, x,y ^ 



{PrfiyW < PTnx)eM^C{T,K,H)\y - x\% 

p-1 



B{\~H,\~H)Y (l+KT) 



whereC{T,K,H) - y-^^rzj^j t^h^h^h) ' 

Proof. By (|4.ip and the Young inequality |6], we have, for all (5 > 0, 

\DyPTf{x)\ < SiPrif log f){x) - {PTf){x){logPTf){x)] + PTf{x)[6logEe^''^]. 
Now let /3(s) = 1 + s{p — l),7(s) = x + s{y — x),s G [0, 1], we have 

£ log(PT/(^))^(7(s)) 

_ p(p-l)PT(/^(^)log/^(^))-(PT/^(^))logPT/^W^^^^^^^ P Dy^^Prf^'^ 



(4.3) 



/32(,) 



Pj,fm 



> 



> 



p 



(ft./'"")logFT/'""(7(«))l - |0„-.-Pr/'""l(7(«)) 

/3^(s) 



his)) 



where we use (14. 3p and choose 6 = for the last inequality, note that Nt is corresponding to 
the direction y — x. 

Next we are to estimate Ee^^^. Since Ee^^^ < ( Ee-s^^^^M^, we turn to the term {Nt). 



(Nt) 



(r(i - H)T)^ Jo 



„2H-l 



{s — r) 



{l+d2bir,Xr)iT -r))dr 



{y — x)'^ds 
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- \Ta-H)T 







dr 



(s — r)2+^ 



' B{l-H,\-H)y {l + KTf 



{y - xfds 



~ \ V{\-H) j T^^2{l-H) 
--: C{T,K,H){y-xf. 



{y - xf 



Therefore, we deduce that 



-I log(PT/^(^))^(7(s)) 



> ^C{T,K,H){y-xf. 

p-1 

Integrating on the interval [0, 1] with respect to s, we get the desired resuh. 

Remark 4.5 To our knowledge, for the results of theorem I3.il and corollary \4.4\ we can not 
decide which is better. 
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